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Abstract 

Upon using the hydrodynamic analog we can derive some families of stationary Beltrami field- 
like solutions from the free Maxwell equations in vacuum. These stationary electromagnetic fields 
are helical and/or column-like once they are represented in a suitable frame of reference. Possible 
dendritic and jelly- fish- like patterns of sprites are demonstrated. 

PACS numbers: 92.60.Pw, 02.40.Yy, 03.50-z, 52.80.Mg, 45.10.Na, 47.32. Cc 
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The majority of the great discoveries in physics have been made by chance (electromag- 
netism, radioactivity, etc.) and only in some rare cases was the experiment preceded by 
theory (lasers, thermonuclear reactors). It seems likely that the problem of sprites belongs 
to the latter one. Since their discovery by Franz et. al. [1], sprites have now been observed 
over thunderstorms all over the world [2-4]. A spectacular manifestation of these sprites 
: the intense transient quasistatic electric (QE) fields of up to ~ 1 kV -m"^, which for 
positive (Cloud-to- Ground) CG discharges is directed downwards, can avalanche accelerate 
upward-driven runaway MeV electron beams, producing brief (~ 1 ms) flashes of gamma 
radiation. These intense fields are large luminous discharges in the altitude range of ~ 40 
km to 90 km, which are produced by the heating of ambient electrons for a few to tens of 
milliseconds following intense lightning flashes. 

Thousands of positive ground flash-sprite associations have been identifled through compar- 
isons with video imaging/optical sensor veriflcation of sprites. Beginning with suggestions 
by Wilson [5], the electrostatic held change of the hghtning flashwas sufficient to exceed 
the dielectric strength of the mesosphere and initiate the sprite. This Wilson mechanism 
for sprites initiated by conventional dielectric breakdown is polarity independent-positive 
and negative changes in charge moment change in excess of the threshold should be equally 
effective in the initiation of sprites. Sprites require a lowly or nonionized medium; they can 
propagate from the ionosphere downwards or from some lower base upwards or they can 
emerge at some immediate height and propagate upwards and downwards [6-7]. 
The first columniform sprites (c-sprites) reported [8] were vertical columns of light that ex- 
tended from about 76 to 87 km and probably less than 1 km in diameter. According to 
Wescott et al. [8] the uncertainty in height determination for their c-sprite measurements 
was approximately 1 pixel or less than 1 km. The standard deviation of the heights for 
the 28 triangulated images was 1.9 km for the top and 1.4 km for the bottom. Additional 
observations of c-sprites, some of which may extend to lower altitudes, have been recently 
reported [9-10]. The bodies of c-sprites usually appear as nearly straight lines, in some cases 
bent slightly, and are composed of bright beads and dark regions. 

Meanwhile, some researchers [7-8] described the phenomenon, called red Sprites^ as a lumi- 
nous column that stretches between 50 and 90 km, with peak luminosity in the vicinity of 
70-80 km. The hashes have an average lifetime of a few milliseconds and an optical intensity 
of about 100 kR. Red Sprites are associated with the presence of massive thunderstorm 
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clouds, and among the most puzzling aspects of the observations is the presence of spatial 
structure in the emissions reported by Winckler et al. [7] . Vertical striations, denominated 
'tendrils', with horizontal size of 1 km or smaller, often limited by the instrumental resolu- 
tion, are apparent in the red Sprite emissions [9-10]. 

As one example, the key features of the previous model [9] are: (1) an ambient conductivity 
profile that falls between a measured nighttime and a measured daytime conductivity; (2) 
an aerosol reduced conductivity in a trail from a meteor that passed through some time 
during the evening, and (3) a cloud-to-ground (hereafter CG) hghtning stroke, with suffi- 
cient charge transfer, subsequent to and occurring within an hour of the development of the 
reduced conductivity trail. This model predicts a temporally brief column of light resulting 
from the conventional breakdown of air in a strong electric field in the observed altitude 
range. The results we shall describe in this short paper might provide an explanation for 
the observed c-sprite phenomena. 
The free Maxwell equations (FME) are 



c dt ' c dt ' 

diEi = 0, diBi = 0, 

where di = d/dxi, Xi — {x,y,z), Cijk is an antisymmetric tensor or symbol. We presume 
that the solutions of FME are 

E{Y,t) = e{r)r{t), and S(r, t) = b(r)^?(i). (1) 

In order to obtain solutions of this system with three constants only and to obtain sinusoidal 
solutions, we also presume that 

-dg{t)/dt _ dJ^{t)/dt _ 

Ht) ~ Q{t) ~ - ^ ^ 

Therefore, the general solutions of above system [11] are 

jr(i) = Co sm{flt - 0), and g{t) ^ Co cos(Oi - 0), (3) 

where Cq and (p are arbitrary constants. Now, equations for e and b become 

Q n 

V X e = -b, V X b = -e. (4) 

c c 
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Above equations give 

V X ^ = — ^, with A^e + h, (5) 

c 

where e is a polar vector (hnked to the electric field E in space, E = eJF(t)), b is an axial 
vector (linked to the magnetic field B in space, B = hQ{t)). 

Now, from the vector analysis and the hydrodynamic analogy (e.g., for a velocity field u, we 
have the vorticity field u; which is defined as V x u) , we know that 

{V xA)\\A or (V xA)xA^ 0, (6) 

here, || means 'parallel to'. In fact, people termed L = (V x ^) x ^ as a Lamb vector [12-13]. 
In generalized fiuid mechanics, we call the fiow which satisfies 

a; X u = (7) 

a Beltrami fiow [11-14], sometimes a screw motion or hehcal fiow [11-13] (here, cv 0). 
With regard to these, vector fields satisfying Eq. (7) are helical or screw fields. It means A 
or the combination of e and b is a helical or screw field. One class of these solutions can be 
obtained and expressed in cylindrical polar coordinates (r, 9,z) as 

u = Aorcos{Q, v = fir + Aqt sm{() , w = -— sin(C), (8) 

where ( = k{z — cot) = kz — 2Qt. When u = {u, v,w) = A is expressed in a frame of reference 
rotating with the mean angular velocity Q {v = Aor sm{()), we have V x u = u; = —ku or 
a; X u = 0. The stationary form of this solution could be similar to e + b (i.e., A). We remind 
the readers that eqns. (8) describe a field that grows proportionally to the distance from the 
cylindrical axis. Therefore, in unbounded vacuum, as there is no dissipation, Agr ~ 0{1) 
as r — > oo should be presumed otherwise they will give some non-localized electromagnetic 
field configurations with infinite energy. Under this consideration, once r — > oo, w is rather 
small (as ^ if k is finite) and this (axial) component (related to magnetic field; will 
be illustrated below) behaves like a weak 'thread'! On the other hand, since Agr ~ 0(1), 
once r ^ 0, there might be a singularity at r = {Aq has been restricted for very larger 
r) . This singularity resembles that of a free vortex (vortex core or the filament of a vortex 
tube). Note that, once the amplitude being replaced by A{t) (with dA/dt — —vAk"^, v 
plays the role of damping, the solution of above equation easily reads A{t) — A(0)e~'^'^^*), 
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we then have u — A{t)r cos{C,) ^ v — A{t)r sin{Q , w — —2A{t)sm.{C,)/k. To demonstrate 
our solutions, we plot those stationary ones of them into following figures. We shall firstly 
present the electromagnetic field in vacuum. Parameters are selected as Aq = 0.45, k — 2.5; 
the maximum spatial range : 0<r<l,0<2;< 27r. The vector field is represented in 
a rotating frame of reference (Fig. 1) with the mean angular velocity Cl = 5.0. We can 
observe the helical and the dendritic features of our solutions and the possible pattern of a 
sprite-hke electromagnetic field. 

The outer and inner isosurface plots (for the total modulus of A and the modulus of uir+vie, 
respectively) are presented in Fig. 1. To make sure there are also possible sprite-like fields 
in our results as there are so many parameters needed for tuning and the 3D vector view is 
sometimes not comprehensive, thus wc only show those isosurface plots below. 
With intensive searching, we can find out, at least, for fl — 5, Aq — 0.45, k = 1.5, there 
exists also possible column-like together with hehcal field for sprites. This pattern is shown 
in Fig. 2 for the same frame of reference ( in a rotating frame of reference). The outer 
(green) surface is for the total modulus of A which is already in a hehcal shape. The inner 
(yellow) surface is for the modulus of uir + vie which demonstrates the possible column-like 
channel. As parameters are tuned into Q = 0.2, Aq = 1.5, k = 2.5, the dendritic features 
of sprites are illustrated clearly in Fig. 3. We can now easily capture the characteristics of 
hehcal or possible helix-like fields in Figs. 1, 2 and 3. These results confirm our present 
approach (cf. images or figures in [8-9]). 
Note that, once we introduce 

e(r) = ^[^(r) - ^(-r)], b(r) = ^[^(r) + ^(-r)], (9) 

considering the definition of Poynting vector S = c(E x B)/(47r), we also can obtain the rel- 
evant results about S (in fact, b(r) oc w as evidenced from above equation). The calculation 
is now only related to ^(r) x ^(— r) instead of e x b. Meanwhile, as the energy density W 
is defined as (E'^ + B'^)/8tt, we can directly interpret W from the magnitude or modulus of 
A from Eq. (8). Meanwhile, from Eqs. (8) and (9), we can obtain the explicit expression 
of the electric (e(r)) and magnetic (b(r))fields. We might have, the electromagnetic energy 
{\S\ oc a-fRo^, a = -{flR^/c) cos(Oi?o/c) + sin(Oi?o/c) and 7 = a - {fl'^Rl/c^) cos(Oi?o/c)) 
within spheres of radius Rq, which is the solution of tan(Qi?o/c) = ^Rq/c, and it is inde- 
pendent on the time. 
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However, once we selected parameters as Q = 5, — 1-2, A; = 0.8, a single sprite appears 
(cf. Fig. 4 in [8] or Fig. 4a in [9] or Figs. 5 and 6 in [15]) which is shown in Fig. 4. The 
other interesting case is for Q = 0.6, Aq = 0.6, k = 0.15 as shown in Fig. 5 : a jelly-fish-hke 
pattern with the downward tendrils-like feature appears. 

To be specific, we like to make a final remark about the characteristic of our Beltrami fields 
: Eq. (8) especially when k is rather small (solutions for /c = violate the Beltrami con- 
dition!). Once A; — > 0, considering the limit of the ty-component in the stationary case, 
w oc z with z being finite. Under this situation, our solutions resemble those results for 
high-electric-field cases [15]. 

In a brief summary, the author likes to point out, the helical and column-like fields illustrated 
above are closely linked to the formation of sprites, however, there will be other unknown, 
say, relativistic, effects during the final stage in the formation of lightning. We shall consider 
other more complicated problems [16-21] in our future works. 
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Q=5,k=2.5,Ao=0.45 
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Fig. 1 Isosurface plots of the solution (e + b) in a rotating frame of reference. 
Possible helical and column-shape electromagnetic fields. Q = 5, Aq = 0.45, k = 2.5. 
(the outer or green surface is for \A\ and the inner or Yghow surface is for [^1^ + viel) 
This result resembles qualitatively that (cf. Fig. 4 in 15]) by Mironychev and Babich 
and those (say, cf. Fig. 1 in jS]) by Wescott et al. It has the dendritic type. 
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Fig. 2 Isosurface plots of the solution (e + b) in a rotating frame of reference. 
Possible column (sprite)-like and helical electromagnetic field. Cl = 5, Aq = 0.45, k 
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Fig. 3 Isosurface plots of the solution (e + b) in a rotating frame of reference. 
Possible sprite- and column-like electromagnetic field. Cl = 0.2, = 1.5, A; = 2.5. 
(the spherical or green surface is for |^|) 
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Q=5,k=0.80,An=1-2 
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Fig. 4 Isosurface plots of the solution (e + b) in a rotating frame of reference. 
Possible sprite and column-like electromagnetic fields. O = 5.0,^40 = 1.2, A; = 0.8. 
The column-like (yellow) isosurface is for e(r). This result resembles qualitatively that 
(cf. Figs. 5 and 6 in by Mironychev and Babich. 




Fig. 5 Isosurface plots of the solution (e + b) in a rotating frame of 
reference. Possible sprite-like electromagnetic field Q = 0.6, Aq = 0.6, k = 0.15. 
Possible jelly- fish- like pattern with the downward tendrils-like feature appears. 



